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ABSTRACT 

We describe high resolution Smoothed Particle Hydrodynamics (SPH) simulations of three approx¬ 
imately M* field galaxies starting from ACDM initial conditions. The simulations are made inten¬ 
tionally simple, and include photoionization, cooling of the intergalactic medium, and star formation 
but not feedback from AGN or supernovae. All of the galaxies undergo an initial burst of star for¬ 
mation at z ~ 5, accompanied by the formation of a bubble of heated gas. Two out of three galaxies 
show early-type properties at present whereas only one of them experienced a major merger. Heating 
from shocks and -PdV work dominates over cooling so that for most of the gas the temperature is 
an increasing function of time. By z « la significant fraction of the final stellar mass is in place 
and the spectral energy distribution resembles those of observed massive red galaxies. The galaxies 
have grown from z = 1 —> 0 on average by 25% in mass and in size by gas poor (dry) stellar mergers. 
By the present day, the simulated galaxies are old (« 10 Gyrs), kinematically hot stellar systems 
surrounded by hot gaseous haloes. Stars dominate the mass of the galaxies up to ~ 4 effective radii 
(« 10 kpc). Kinematic and most photometric properties are in good agreement with those of observed 
elliptical galaxies. The galaxy with a major merger develops a counter-rotating core. Our simulations 
show that realistic intermediate mass giant elliptical galaxies with plausible formation histories can 
be formed from ACDM initial conditions even without requiring recent major mergers or feedback 
from supernovae or AGN. 

Subject headings: galaxy formation: general — galaxy formation: elliptical — methods: numerical 


1. INTRODUCTION 

There are many puzzles encountered in understand¬ 
ing the formation and evolution of elliptical galaxies and 
the spheroidal components of spiral galaxies. On the 
one hand a naive reading of the hierarchical theory of 
structure formation in a ACDM universe would argue 
that, since massive halos form later than less massive 
ones, massive ellipticals which reside in the centers of 
these massive halos should also form late. But there 
is strong observational evidence that old, massive, red 
and metal rich proto-ellipticals are already in place at 
z = 2 — 3 and that present day early-type galaxies 
formed most of their stars well before a redshift z = 1 
(ISearle e t al.lll973tlBrinc hmarm fc Ellisll2000b iTreu et all 
120051: Ivan der Wei et, all [20051 1. We also know from ev¬ 
idence dating to the 1970’s that current rates of star 
formation in these systems are quite low , the rates i n¬ 
creasing sharply into the past as z 1 (Butcher fc Oemlerl 
1 197 8l fDressledll980l : lltakos fc Schombertl l 19951. Further, 
with regard to the evolution of this population, there 
appears to be a significant increase in the total stellar 
mass in these old objects from z = 1 to the present 

dBell et alJl2004lprorv et abll2004 IConselice et al.ll2005l : 

iFaber et al.l 1200511 . but this cannot easily be accounted 

for by the fading o f younger, bluer s t ar forming galaxie s 
present at z = 1 (|Drorv et al.ll2004 IFaber et al l 120051) . 
This has led to the plausible and popular idea that ‘dry 
merging’ (i.e. merging of predominantly stellar systems) 
among elliptical system s (ITran et al.ll2005l Ivan Dokkuml 
1 2 005t [Bell et al.l f2006bl fah pushes more and more of the 
mass over the observational cutoff to provide an increased 
number of 0.5-2L* galaxies. 


This attractive idea gives reasonable explanations 
for some observed correlations of the mass of el¬ 
lipticals and their kinematics, sizes or isophotal 
shape dKhochfar fe Burkertl 120051 : iNaab et ahl l2006bl 
iBovlan-Kolchin et al.ll200fil) . However, noting that there 
are tight relations among luminosity, color, age and met- 
alicity, puts severe limits on the amount of late dry merg¬ 
ing that could acceptably occur without destroying the 
aforementioned constraints. 

On the theoretical side it is now possible to perform 
accurate high resolution simulations of the gravitational 
evolu t ion of the dark matt er distribution (jMoore et al.l 
Il998t ISnringel et al.l l2005bt ). In contrast, the numeri¬ 
cal simulation of galaxy formation, including a hydrody¬ 
namic treatment of the baryonic component, is still in 
its infancy. Very few high resolution simulations from 
realistic cosmological initial conditions have been done 
so far and most of these have concentrated on the for¬ 
mation o f disc galaxies rather than early-type spheroidal 
systems (ISommer-Larsen et al.l l2003t lAbadi et al.l 120031 

iGovernato et al.ll2004 iRobertson et al.ll2004lh 

This is surprising since spheroidal systems are of in¬ 
terest in their own right, as they contain more than 
half of the total ste llar mass in the local universe 
(|Fukugita et al.l 1199811 . The most massive galaxies 
known, the giant ellipticals, are spheroidal systems w hich 
predom inantly consist of old stars (see e.g. iThomas et alJ 
l2005all and so must have formed at high redshift. They 
are therefore likely to be good probes of galaxy assem¬ 
bly, star formation and metal enrichment in the early 
Universe. 

Most numerical work on early-type galaxy forma- 
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Fig. 1.— Circular velocity curves for galaxy A at four different 
numerical resolutions: 40 3 , 50 3 , 100 3 , and 200 3 SPH particles and 
collisionless dark matter particles, respectively. Note how the ro¬ 
tation curves become increasingly flat as the resolution increases. 
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Fig. 2.— Star formation rate (SFR) histories, computed from 
stellar ages, of galaxy A versus lookback time at four different 
numerical resolutions: 40 3 , 50 3 , 100 3 , and 200 3 SPH particles and 
collisionless dark matter particles, respectively. There is a strong 
trend that the low redshift star formation rate is reduced in higher 
resolution simulations. 

tion has used ei ther very idealized initial conditions 
dKobavashil 120041) . or had insufficient spatial and mass 
resolution to resolve the internal structur es of galaxies 
(jKawata fc Gibsoiil 120031 ISaiz et all l2fl04H An excep¬ 
tion is the simulation discussed bv iMeza et all (2003). 
These authors used an SPH simulation that included 
feedback from supernovae to follow the formation of a 
single spheroidal galaxy self-consistently from CDM ini¬ 
tial conditions. The spatial resolution of this simulation 
was high enough to resolve the region within an effective 
radius for a typical real elliptical galaxy. However, the 
final stellar system formed in this simulation was far too 
dense, with an effective radius about an order of mag¬ 
nitude smal ler than rea l ellipt ical galaxies of the same 
brightness. IMeza et 'ail 1 20033 speculate that this dis¬ 
crepancy may be a consequence of their star formation 
and feedback algorithm and that it might be possible to 
produce less concentrated systems if more aggressive stel¬ 
lar feedback were implemented to prevent star formation 
in high density sub-units at high redshift. 


In this paper we present high resolution hydrodynam- 
ical simulations based on cosmological initial conditions 
admitting, intentionally, only bare-bones prescriptions 
for the physics involved (e.g. no “feedback” from su¬ 
pernovae or AGN), to see if some resolution of these 
paradoxes can be derived. The rationale for simplyf- 
ing the simulations is straightfoward. As we will show 
in Section [3l at least 100 3 gas particles are required for 
hydrodynamical simulations to ‘converge’, making them 
very expensive to run. With present generation comput¬ 
ers it is impossible to run a large ensemble of simulations 
of this size to properly explore a huge parameter space. 
Our point of view, therefore, is to keep the physics of the 
simulations as simple as possible and to get an under¬ 
standing of the behaviour of a simplified problem before 
investigating additional complexities such as supernova 
and AGN feedback. Our goalin this paper, therefore, 
is to investigate the formation of a number of massive 
(~ M*) isolated galaxies starting from realistic initial 
conditions and to see how variable the final systems and 
whether they resemble real galaxies. 

What we find can be summarized simply. The ini¬ 
tial cooling and collection of cold gas from infalling 
smaller scale perturbations happens very rapidly and eas¬ 
ily within the most massive systems leading to a very 
rapid burst of star formation beginning at z ~ 6 and 
then falling off exponentially on a time-scale of roughly 
1.5Gyrs. This phase is terminated as the star forming re¬ 
gion is enveloped in an expanding hot bubble which pre¬ 
vents new, infalling cold gas mass elements from reach¬ 
ing the central regions. This early phase - reminiscent in 
some respects to a modernized version of the “monolithic 
collapse” picture of gal axy formation, prominent in th e 
1960’s through 1980’s (|Larsonl[l974 Ivan Albadal 1 1982th 
produces a sequence of objects with effective radii of 1- 
2kpc, which might satisfy the tight relations observed 
among the red metal rich old cores of elliptical galaxies. 

For simulations with early-type properties at z = 0 
stellar accretion or mergers (the choice of the appropri¬ 
ate term for the process is arbitrary and dependent on 
the relative mass of the infalling stellar objects) add to 
the growing stellar envelope of relatively blue, old and 
metal poor stars. This accounts for the growth in size 
seen for ellipticals as well as the growth in total mass in 
the time frame z = 1 —> 0. Furthermore, since this as¬ 
sembly of added stellar mass is not accompanied by much 
in situ star formation (the coterminously added hot and 
cold gas being simply added to the expanding hot bubble 
surrounding the central galaxy) the growth of the stellar 
population occurs without the presence of young stars. 
Existing work has shown that the tightness of the ellip¬ 
tical color-magnitude r elation puts strong constraints on 
dry m erger scenarios (1 Bower et al.1 Il998t iKavirai et al.l 
f2005h . But we find that in our simulations that minor 
mergers or accretion events do not typically add much 
stellar mass to the central region, because of the rela¬ 
tively large angular momentum of the infalling systems - 
and thus the central, tight color-magnitude relations are 
not expected to be overly perturbed even if as much as 
40% of the total final mass is added during this phase. 
As an additional byproduct of this assembly scenario, we 
find it readily understandable how so much of the stellar 
mass in these systems resides in regions of such low den¬ 
sity, where star formation would always have been diffi- 
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cult to contemplate. The stars seen at R ~ R e were not 
born there, but rather were formed in the central regions 
of much lower mass systems that have been accreted and 
shredded in the envelopes of the giant elliptical galaxies. 

2. INITIAL CONDITIONS AND SIMULATIONS 

The initial conditions of the ACDM simulation as¬ 
sumed scale-invariant adiabatic fluctuations. The post¬ 
recombination power spectrum for the C DM cosmology 
was b ased on the parameterization of Efstathiou_et - ^L 
(1992) with T=0.2. We use a WMAP1 (|Spergel et al. 
I2003fl cosmology with a slightly lower Hubble parameter 
of h = 0.65 (= H o =100h kms~ 1 Mpc” 1 ) with tjg=0.86, 
fb = Hf,/fl m =0.2, fl 0 =0.3, and A 0 =0.7. 

In a low resolution, dark matter simulation (128 3 dark 
matter particles with a box size of Lh nx = 5 0 Mp c 
using the (AP 3 M) N-body code of ICouchmanl (|1991l) ). 
We selected virialized halos with masses in the range 
7 x 10 n M 0 < Mhaio < 3 x 10 12 Mg in a low density en¬ 
vironment such that the nearest halo with Mhaio > 2x 
10 u Mq is more than one Mpc /h away. To re-simulate 
the target halos at high resolution we increased the parti¬ 
cle number to 100 3 gas and dark matter particles within 
a cubic volume at redshift z = 24 containing all particles 
that end up within the virialized region (conservatively 
we assumed a fixed radius of 0.5 Mpc) of the halos at 
z = 0. From this sample we have selected three systems 
with final stellar masses of ss M* for elliptical galaxies 
in the range of 9.6 x 1Q 10 M(d < M < 1.6 x IO^AMq 
( C ole et all 200 1) iKochanek et all[2001 : [Bell et afll200 3fi 
for the cosmology chosen here. For one halo (halo A) the 
resolution was increased to 200 3 particles. Additional 
short wavelength perturbations were included to account 
for the missing small-scale power below the Nyquist fre¬ 
quency of the low resolution simulation. The tidal forces 
from particles outside the high resolution cube were ap¬ 
proximated by increasingly massive dark matter particles 
in 5 nested layers of lower and lower resolution. None 
of the galaxies was contaminated by boundary particles 
within their virial radii. 

The simulations were run with GADGET-2 (jSpringell 
120051) on COSMOS, a shared-memory Altix 3700 with 
1.3-GHz Itanium2 processors hosted at the Department 
of Applied Mathematics and Theoretical Physics (Cam¬ 
bridge) and on two similar machines at Princeton Univer¬ 
sity and the University Observatory in Munich. We used 
a fixed comoving softening until z = 9, and after this the 
softening remained fixed at the same values in physical 
coordinates. The simulation parameters for each halo are 
given in Tab[3l The amount of stars formed for each gas 
particle is proportional to the mass of the star forming 
gas component and inversely proportional to the star for¬ 
mation timescale t*. This timescale is i* = t q yjp C iit/p, 
where t\ j is the characteristic star fo rmation timescale 
which we set to = 1.5 Gyr dSprineel fc Hernauistl 
120031) . Within a factor of a few the exact value of tj) 
does not change the results. 

The density threshold for the onset of star forma¬ 
tion was set to p rr it = 7 x 10 -26 gcm~ 3 following 
iNavarro fe White! (1 19931 ). In addition we require an over¬ 
density contrast of A > 55.7 for the onset of star forma¬ 
tion in order to avoid spurious star formation at high red- 
shift. We included an uniform UV backgroun d radiation 
field of a modified lHaardt fe Madaul (j 19961 ) spectrum, 



Radius [kpc] 

Fig. 3. — Time evolution of the gas temperature profile from 
z = 5 to z = 0 (from bottom to top) for halo A (200 3 resolution). 
The average temperature of the gas is steadily increasing. At the 
end of its initial formation phase &t z & 2 the galaxy is surrounded 
by a halo of hot gas heated to the virial temperature. 



Fig. 4. — Spectral energy distribution (SED) of the galaxy A 
(200 3 ) at 2 = 1 (black line) with solar metalicity for all stars. 
We have assumed that stars younger than 10 7 yrs are obscured. 
For comparison we show the SED for a 3.5 Gyrs old simple stellar 
population (SSP) with so lar metalicity (grey dotted line) using the 
IBruzual Sz Chariot! (120031) models with a Salpeter IMF. 

where reionization takes place at 2 : ~ 6 dDave et al.lll999D 
and the intensity o f the UV b ackgro und field peaks at 
z ~ 2 — 3. Followine lKatz et~akl (119961) the radiative cool¬ 
ing and heating rates were computed assuming that the 
gas is optically thin and in ionization equilibrium. The 
self-s hielding of star-forming clouds from the UV flux 
(e.g. ISusa &; Umemural 120001 iZheng fe Miralda-Escudel 
I2002D was not included in the simulations. The abun¬ 
dances of the different ionic hydrogen and helium species 
were then computed by solving the network of equilib¬ 
rium equations self-consistently for a specified value of 
the UV background radiation field. The late star forma¬ 
tion history is relatively insensitive to minor changes in 
the UV background field, whereas an increase in the UV 
background at high redshift might soften the star forma¬ 
tion peak. A detailed investigation of the effect of the 
UV background on the simulations presented here will 
be discussed in a future paper. 

3. A MASSIVE GALAXY AT Z=1 
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TABLE 1 

Properties of galaxy A at r < r vir 


Resolution 

M vir 

A^stars 

A'/gas 

A^dark 

r . (b) 

' vir 

Umax 

TTistars 

^dark 

Cstars 

e dark 

40 3 

242 

25.3 

20.2 

197 

434 

341 

161 

1288 

0.625 

1.3 

50 3 

241 

27.0 

20.1 

193 

433 

356 

82.5 

659 

0.5 

1.0 

100 3 

230 

23.8 

20.9 

184 

427 

270 

10.5 

82 

0.25 

0.5 

200 3 

225 

24.9 

18.4 

182 

424 

232 

1.3 

10.3 

0.125 

0.25 


Note. — (a) Total masses M inlO lo M 0 ; (b) Virial radius in 
10 5 A/g; (e) Gravitational softening lengths in kpc 


TABLE 2 

Properties of galaxy A at r < 30kpc 


Resolution 

A^stars ^ 

.Mgas 

A^dark 

40 3 

16.3 

1.0 

27.7 

50 3 

16.6 

1.0 

28.7 

100 3 

12.2 

0.6 

20.4 

200 3 

11.8 

0.5 

23.1 


Note. — (a) Total masses M inlO lo M 0 

To investigate the effect of numerical resolution on the 
cosmological simulation of an individual galaxy we have 
simulated our primary example, halo A, at four different 
resolutions using 40 3 , 50 3 , 100 3 and 200 3 SPH particles 
and the same number of dark matter particles, respec¬ 
tively. The simulation parameters are given in Tab. [I| 

The virial radius and the enclosed total mass in 
baryons and dark matter of the galaxy changes weakly 
with increasing resolution (see Tab. [I]). However, the 
amount of baryons concentrated to the central 30 kpc at 
the higher resolutions is reduced by more than 20% com¬ 
pared to the lowest resolution (see Tab. [2]). In addition 
the circular velocity profiles v% = GM[r)/r become flat¬ 
ter with increasing resolution and the peaks are shifted 
to smaller radii (Fig. [[]). 

The detailed star formation history of the galaxies also 
changes with numerical resolution. In Fig. [2]we compare 
the SFR histories of the 40 3 , 50 3 , 100 3 , and 200 3 simula¬ 
tions. The numerical resolution has a weak influence on 
the SFR at redshifts higher than z = 1. In all cases the 
SFRs peak at z « 5 at rates of ~ 6OM 0 /yr. However, 
at the 200 3 resolution the peak of the SFR is shifted to 
higher redshifts. At redshifts below z = 1 the low reso¬ 
lution galaxies, e.g. 40 3 and 50 3 particles, have signifi¬ 
cantly higher star formation rates (and gas fractions, see 
Tab. [2]) than the 100 3 and the 200 3 simulation which has 
the lowest present day star formation rate. The higher 
star formation rates (and larger final baryon fractions) 
of the low resolution simulations are caused by a contin¬ 
uing infa ll of cold gas, even a t lo w redshifts. It has been 
shown bv lAgertz et al.l (|2006l i and I Junk et al.l (|2006l i that 
SPH, especially at low resolution, has a limited abillity 
to resolve Kelvin-Helmholtz instabilities: Instead of be¬ 
ing destroyed a cold gas cloud moving in a hot gaseous 
halo can be artificially stabilized and sink all the way to 
the center. In our simulations this artifical effect is sup¬ 
pressed at higher resolution as the infalling clumpy gas 
is better resolved and can be stripped and dispersed in 
the hot gaseous halo of the host galaxy. 

A comparison simulation at the 100 3 resolution with 
a longer star formation time-scale of ijj = 4.5 Gyr is 
similar to the = 1.5 Gyr simulation at high redshifts 
but shows significantly less star formation at low redshift. 
These results implicate that the resolution has a strong 


(c) Maximum circular velocity in km/s; (d) Particle masses m in 

effect and the actual choice of the time-scales has a weak 
effect on the galaxy properties. As the present day color 
of the galaxy is very sensitive to recent star formation it 
is very problematic to use it as a diagnostic tool without 
considering other galaxy properties. 

We have used the highest resolution (200 3 ) simula¬ 
tion to study the star formation in some more detail. 
The stars form in clumps of cold gas that have col¬ 
lapsed on small scales. As the cold gas is depleted, 
the SFR declines rapidly. The diffuse gas in the form¬ 
ing halo, which does not reside in cold sub-clumps, is 
heated by shocks and forms an expanding bubble of hot 
gas which reaches virial temperature at 2 < z < 3. 
As an example we show the evolution of the gas tem¬ 
perature of halo A in Fig. [3] Although during this 
phase the cooling time of the hot gas is shorter than 
the Hubble ti me, the heating time of the infalling gas is 
even shorter (| Ciotti fc Ostri kcr 2001 : iBirnboim fc Dekell 
120031 [Sazonov et al.l 120051 : iDekel fc Birnboimll2006fb Af¬ 
ter z = 2 there is an increasing contribution to heating 
from -PdV work and entropy addition via weak shocks 
which becomes dominant towards z = 0. In this paper we 
focus on the assembly history and the observable proper¬ 
ties of the simulated systems. A full quantitative analysis 
of the heating processes is beyond the scope of this pa¬ 
per and will be presented in Johansson et al. (2006, in 
preparation). 

Galaxy A (200 3 ) has assembled 80% of its final stel¬ 
lar mass by z = 1 (see Section [1. At this time the 
stellar mass within 30 kpc (we use this fiducial value as 
the stars in ellipticals typically can not be observed at 
larger radii) is M* = 9.8 x 1O 1O M 0 and the mean age of 
the stellar population is 3.9 Gyrs. We constructed the 
spectral energy distribution (SED) of the galaxy at z = 1 
using the spectro-photom etric stellar evolution models of 
(|Bruzual fc Charloo 1200 3) assuming a mean stellar met- 
alicity close to solar Z = Zq = 0.02 and a Salpeter IMF 
(Fig. SI)- Qualitatively the SED does not change for 
small changes in metalicity. We also assumed that re¬ 
cently formed stars younger than 10' yrs are obscured 
and do not contribute to the total flux. No further cor¬ 
rections for reddening, e.g. for absorption by dust, have 
been applied. The SED looks remarkably sim ilar to rest- 
frame SEDs of ERO’s (|Moustakas et al.1l200 ll) . However, 
in the absence of reddening corrections, the redshift cor¬ 
rected observed colors are not quite extreme enough at 
R — K = 4.6 and I — K = 3.2. With an observed magni¬ 
tude of mu = 19.3 and the above red colors the object 
could al ternatively be classified as a Very Red Objec t, 
see e.g. (lMcCarthvll2004t IVaisanen fe Johanssonll2004jl . 

4. ASSEMBLY HISTORIES FOR THREE GALAXIES 

The similarity in the star formation and assembly his¬ 
tory as well as in the present day properties of the 100 3 
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Fig. 5. — Top: Dark matter mass accretion history of the three 
(100 3 ) halos within their virial radii. Halo E accretes smoothly 
after z = 1. Bottom : Evolution of the halo spin parameter A . The 
merger events for halos A and C at z 0.7 can be clearly identified 
by the sudden increase of the spin. 

and 200 3 simulations for galaxy A encouraged us to be¬ 
lieve that at a resolution of 100 3 particles (and the cor¬ 
responding spatial and mass resolution) the simulations 
are numerically well enough resolved to run more halos 
and find out whether the above results are typical. With 
current hardware and software it is difficult to run many 
halos at 200 3 or even higher resolution. 

The dark matter mass accretion history and the spin 
evolution of the three halos is shown in Fig. [5j We 
computed the mass within the virial radius defined by 
an enclosed mean overdensity of 200. By z = 1 halos A 
and C have gained « 40% of their final mass. Thereafter 
the halos continue to accrete dark matter and halo A has 
a merger with a mass ratio of 6.5:1 at z « 0.6, halo C 
has a major merger with a mass ratio of 3.5:1 at z ~ 0.8. 
Halo E has already assembled more than 60% of its mass 
by z = 1 and undergoes an almost equal-mass merger at 
z = 1.6. Thereafter the halos grow until z = 0 without 
any indication for a significant merger event. No halo 
experienced a major merger within the last 4 Gyrs. 

The evolution of the halo spin parameter A defined as 


V2M vir v c Tvir 

where v c is the circular veloc ity at the virial ra dius 
rvir and M v i r is the virial mass (Bullock ct al. 200T1). is 
shown in the bottom plot of Fig. [5j Halo E had gained 
significant spin during the z = 1.6 major merger event 
which can clearly be identified by the sudden jump in 
A . This merger indu ced increase of spin is in agreement 
with the findings of lVitvitska et al.l (|2002f ) based on dark 
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Fig. 6.— Top: Baryonic mass accretion history within the virial 
radius of the galaxies subdivided in stars and gas. Bottom: Star 
formation rate histories computed from stellar ages at z = 0. The 
dashed line indicates an evolution ex exp (—(t — to)/ T ) with to = 
14.5 Gyrs and r = 1.5 Gyrs for comparison. 

matter simulations. By z = 1 the spin of halo E is larger 
than A’ = 0.05 and decreasing smoothly towards A = 
0.03 at z = 0. In contrast, Halos A and C have low 
spin at z = 1 and thereafter gain significant spin during 
the merger events. The final spin values of all halos, 
A* < 0.04, is below the average predicted by dark matter 
simulations, however, their time evolution is consistent 
(|Vitvitska et al.ll2002l) . 

In the top row of Fig. [6] we show the accretion his¬ 
tory of the baryonic matter subdivided into gas and stars 
within the virial radius of each halo. Halo A has assem¬ 
bled most of its stellar mass by z = 1 whereas halos C 
and E grow significantly at lower redshifts. By 2 = 0 
halo A has almost half of the baryons in gas with a ratio 
of stellar mass to gas mass of M stars /M gas = 1.1. Halo 
E, at a similar mass, has significantly more baryons in 
stars with M a t ais /M gas = 1.6. We derive the star for¬ 
mation rate (SFR) histories of the galaxies (bottom plot 
of Fig. E) using the stellar ages at 2 = 0 of all stars 
within 30 kpc (again, we use this fiducial value as the 
stars in ellipticals typically can not be observed at larger 
radii). All galaxies start forming their stars in a burst at 
2 ~ 3 — 5 with peak SFRs of ~ 40 — 6OM 0 /yr and un¬ 
dergo an intense phase of merging (see the random walk 
of the halo spin parameters in Fig. 0 during which the 
SFRs decline exponentially towards 2 = 1. At 2 = 1 the 
galaxies sustain star formation at rates of 3 — 7M 0 /yr. 
For halos A and C the star formation rate drops to zero 
soon after their last merger event. About 2 Gyrs later 
small amounts of the previously heated gas is cooling at 
the center of these galaxies and the SFR is rising again 
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Fig. 7. — Accretion history for galaxies A,C and E of stars that 
have been added onto the central physical 30 kpc (asterisk), stars 
that have formed in the galaxy in situ (open diamond) and the 
sum of the two (filled squares). 


to ss 2M 0 /yr. Halo E, which did not undergo any sig¬ 
nificant merger after z = 1 sustains an almost constant 
star formation rate until z = 0. 

The above star formation histories were derived from 
the present day fossil record of the stars in the galaxies. 
They indicate the time but not the location the stars 
were born. The stellar mass of the galaxies is steadily 
increasing with time both by stars that were born from 
gas in situ in the galaxy as well as stars that have formed 
outside the galaxy and later on have been accreted. In 
Fig. [7] we show the time evolution of stellar mass within 
a fixed physical radius of 30kpc separated into stars that 
were born within the galaxies and stars that have been 
accreted. By z = 1 galaxy A has already assembled 80% 
of its final stellar mass when about 60% of the stars have 
formed inside the galaxy whereas 40% have been accreted 
by mergers. From z = 1 to z = 0 further mass assembly 
is dominated by accretion of stars which are accreted in 
small subunits. The most massive satellite is accreted at 
z ~ 0.2 and, after it has been stripped, has 10% of the 
mass of the host galaxy. It had originally been accreted 
onto the halo in the merger event at z = 0.6 that is 
described above. Galaxy C experiences a major ss 3:1 


merger at z = 0.6. This merger can be identified as 
a gas poor or dry merger as both progenitors were gas 
poor and the merger was not accompanied by a burst 
of star formation (see Fig. [6]). Galaxies A and C show 
on average an increase in mass by ‘dry’ mergers of 24% 
and only a few percent by in situ star formation from 
z = 1 — > 0. The situation is different for galaxy E which 
grows predominantly by in situ star formation from 2 = 
1 —> 0 (see Fig. [9]). 

In Fig. [5] we show the time evolution of the projected 
cylindrical half-mass radius (again within fixed physical 
30kpc) for the galaxies as a whole and for the in situ 
and accreted population separately. The stars formed in 
situ have typical half mass radii of 1-2 kpc which changes 
only weakly with redshift. The accreted stars have sig¬ 
nificantly larger half-mass radii that lead to an increase 
of the total half-mass radius for galaxies A, C and E as 
the accreted mass is increasing towards the present time. 
The situation is different for galaxy E: the in situ half 
mass radius is constant as well, however, the mass in in 
situ stars is increasing. As a result the half mass radius of 
the galaxy is decreasing similar to adiabati c contraction 
(iBhimenthal et al]l!986 : l.lesseit, et, al.ll2002l ). 


5. PRESENT DAY PROPERTIES 

We have summarized the properties inside the virial 
radius of the galaxies in Tab. [3] In Fig. [lU] we 
show the circular velocity curves for all galaxies normal¬ 
ized to their projected stellar half mass radii. These 
seem very rea sonabl e, in marked contrast to the re¬ 
sults of e.g. iMeza et all ( 2003 ). The galaxies are 
dominated by luminous matter inside 10 kpc and the 
fraction of dark matter within the half-mass radii is 
0.18 < /dm ~ 0.38. This is consistent with most find¬ 
ings based on observations and dynamical modeling that 
early-type galaxies ar e dominated by l umin o us matter 
in their inner parts jE fstathio u et al.l fl 982: Rix et all 
1997t iRomanowskv et al.ll200.3l iTreu fc Kooomansll2004 

Thomas et al.ll2005bD . However, dynamical modeling of 
observed ellipticals also indicates that dark matter be¬ 
comes dominant beyond 2-3 half-mass radii whereas for 
the simulated galaxies A,C, and E dark matter becomes 
dominant at a factor of two larger radii. 

The properties of the central stellar component of the 
galaxies are summarized in Tab. |4] At present galaxies A 
and C have very old stellar population with mean ages of 
about 11 Gyrs, galaxy E which had significant recent star 
formation is more than one Gyr younger. The final stellar 
masses within 30kpc are in the range of 10.8 x 10 lo Mg < 
M* < 13.2 x 10 lo M Q . We ha ve computed the luminosi - 
ties of the galaxies using the iBruzual fc Gharlotl (2003:) 
spectro-photometric models assuming a Salpeter IMF 
and solar metalicity. As all galaxies had recent star for¬ 
mation on various levels (see Fig. m, they are shifted 
blue-ward from th e colo r magnitud e relation of cluster 
ellipticals (see e.g. iMcIntosh et al.l 12005). However, for 
galaxies A and C the effect is not very strong and any 
process that would suppress central star formation within 
the last 10 8 yrs would shift the galaxies (with the excep¬ 
tion of galaxy E) to the observed relation. Higher numer¬ 
ical resolution might result in a similar effect, as galaxy 
A at 200 8 resolution falls right on the color-magnitude 
relation. However, it is not unusual for field ellipticals 
(E+A phenomenon) to actually show residual star for- 
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TABLE 3 

Galaxy properties for r < r vir 


Galaxy 

M vil » 

Af stars 

A/gas 

Afdark r v ir (b) 

Umax 

< A‘ > W 

fhot (e) 

in Stars 

^dark 

^stars ^ 

e dark 

A 

229 

23.8 

20.9 

184 

426 

270 

0.022 ± 0.012 

0.99 

10.5 

82.5 

0.25 

0.5 

C 

158 

18.0 

13.2 

127 

377 

257 

0.034 ± 0.012 

0.99 

8.5 

69 

0.25 

0.5 

E 

171 

21.3 

13.1 

137 

387 

307 

0.039 ± 0.010 

0.98 

8.5 

69 

0.25 

0.5 


Note. — (a) Total masses M in 10 10 Mq : (b) Virial radius in kpc; (c) Maximum circular velocity in km/s; (d) Time averaged halo spin 
parameter; (e) Fraction of hot gas (T > 10 s K); (f) Particle masses m in 10 5 Mq; (g) Gravitational softening lengths in kpc 


TABLE 4 

Galaxy properties for r < 30kpc 


Galaxy M sta rs 

Af^as 

hot (b) 

G/2 (c) 

<age>* 

M B (e) 

SB eff W 

Teff (g) 

/DM (h) 

e 0> 

d4 ^ 

a 0 < k > 

v -d) 

b maj 

< h 3 > 1 

» T <r 4 

A 

12.3 

0.64 

0.95 

2.2 

11.5 

-20.18 

20.5 

2.7 

0.38 

0.18 

0.6 

153 

43 

-0.05 

0.08 

C 

10.8 

0.82 

0.90 

2.5 

11.2 

-20.45 

20.3 

2.8 

0.35 

0.21 

1.1 

145 

84 

-0.09 

0.77 

E 

13.2 

0.84 

0.88 

1.4 

9.8 

-21.06 

19.7 

1.9 

0.18 

0.5 

2.5 

180 

225 

-0.05 

0.13 


Note. — (a) Total masses M in ! O 10 A/q : (b) Fraction of hot gas (T > 10 5 K); (c) Projected stellar half-mass radius; (d) Average stellar 
age at z=0; (e) Absolute magnitude in B-band; (f) Effective surface brightness in B-band in mag arcsec -2 ; (g) Fitted effective radius; (h) 
fraction of dark matter within r e g; (i) Ellipticity at r e fj; (j) Isophotal shape parameter; (k) Central stellar velocity dispersion in km/s; (1) 
Major axis velocity in km/s; (m) LOS YD assymetry parameter; (n) Triaxiality parameter 
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Fig. 8.— Time evolution of the projected spherical half-mass 
radius for galaxies A,C and E of all stars (filled squares) and the 
accreted (stars) and in situ stars (open diamonds) within a fixed 
radius of physical 30 kpc. 


Fig. 9. — Accreted stellar mass and stars formed in situ as a 
function of redshift for the three galaxies. From z = 1 —> 0 galaxies 
A and C grow by gas poor mergers whereas galaxy E grows by in 
situ star formation. 
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Fig. 10.— Circular velocity curves for all galaxies versus radius 
normalized to their projected stellar half-mass radius. All galaxies 
are dominated by luminous matter inside two half-mass radii. 




Fig. 11.— K-band surface brightness profiles for galaxies A,C, 
and E. The profiles inside ps 2 softening lengths (shaded area) were 
excluded from the fits as they are likely to be influenced by force 
softening. The error bars are given by the Poissonian error of the 
particles in the fitted regime. The errors for the Sersic index n and 
the effective radius r e have been determined by the bootstrapping 
method. Galaxies A and C can be fitted well by a single Sersic 
function whereas galaxy E can only be fitted by a two component 
fit (see text). 


mation (|Blake et al.ll2004 ISchawinski et al.ll200(il) . 

The simulations had sufficient resolution (a few 10 5 
particles within 30kpc) to analyze the kinematic and pho¬ 
tometric properties in more detail. In Fig. [lT]we show 
the edge-on K-band surface brightness profiles of the 
galaxies and the best fitting Sersic function o r the best 
two c omponent Sersic+exponential fit dNaab &; Truiillol 
120061) . With Sersic indices of n = 4.3 and n = 2 and ef¬ 
fective radii of 2.7 kpc and 2.8 kpc galaxies A and C can 
best be fitted by a single Sersic function and have a very 
reasonable concentration in good agreement with ob¬ 
served intermediate mass ellipticals dTruiillo et al.ll2004ll . 
Galaxy E can not be fitted with a single component. It is 
dominated by a central ’bulge’ component with a Sersic 
index of Ub = 0.5 and a size of r e s = 1.7kpc and an 
outer exponential with a scale length of hd = 4.8fcpc. 
The bulge-to-total ratio is B/T = 0.9. In fact the central 
component is a compact exponential disk seen edge-on. 

We have investigated the LOSVDs of the galaxies 
by binning the velocities on a slit along the appar¬ 
ent long axis of each projected galaxy. Subsequently 
we parameterized deviations from the Gaussian shape 
using Gauss-Hermite basis functions (jGerhardl 119931 
Ivan der Marel fe Franxl 119931) . The kinematic parame¬ 
ters of each profile (vat, er^t, h 3 , / 14 ) were then determined 
by least squares fitting as in iNaab et al. (l2006a l. As an 
example we show the fitted kinematical profiles for the 
edge on projection of galaxies A and C in Figs. [T2land[l3l 
As they are affected by higher order moments we show 
the real line-of-sight velocity and dispersion in addition 
(dashed lines). Both galaxies show significant rotation. 
Galaxy C has a counter-rotating core that has formed 
during the major merger event. The velocity dispersion 
profiles are peaked towards the central region and drop 
in the outer parts. The drop at the very center is caused 
by stars that have formed late from the cooled halo gas. 
A similar behavio r is seen in some field ellipticals like, 
e.g. NGC2768 (|Emsellem et al.ll2004f) and could have the 
same origin. The asymmetry of the LOSVD measured by 
/13 is always anti-correlated to the line-of-sight velocity, 
whereas the LOSVD is more peaked than a Gaussian at 
the center, indicated by a positive value of / 14 . The errors 
have been computed applying the bootstrap algorithm. 
Even for the higher order moments they are reasonably 
small within the effective radius. 

In Fig. [L4] we show the location of 500 random pro¬ 
jections of galaxies A,C, and E in the vm/aat -/13 plane. 
The shaded area indicate the 90% probability to find pro¬ 
jected galaxy with the given properties. All projected 
galaxies show steep leading wings in their LOSVD (anti¬ 
correlated /13 and The location of the galaxies is in 
good agreement with observations of elliptical galaxies 
indicated by the s ymbols in Fig. [14] (|Bender et all 1 19941 : 
INaab et al.ll2006ai ). Galaxy E is the fastest rotator with 
'Dit/CTfit ~ 2 . 

We determined the characteristic ellipticity e e s, the 
isophotal shape a4 e ff (a4 e ff < 0: boxy isophotes; a4 e g > 
0 : disky isophotes) as well as the central velocity dis¬ 
persion do, the maximum rot ation velocity iWi, and the 
triaxiallity parameter T as in lNaab fe Burkern (120031 ). In 
Fig. [I5]we indicate the location of the simulated galaxies 
A,C, and E in the plane defined by the anisotropy pa¬ 
rameter 'e*iaj/<r 0 an d the isophotal shape a4 e ff. We have 
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Radius [r/r e ] 

Fig. 12.— Analysis of the LOSVD of galaxy A as measured along 
a slit aligned with the major axis of the moment of inertia tensor 
of the main stellar body (open squares). The fitted local velocity 
Ufit, local velocity dispersion , v TO t /crfit» # 3 , and H 4 are plotted 
versus radius. The dashed lines show the true line-of-sight velocity 
and dispersion, respectively. The individual error-bars were derived 
by bootstrapping. 

computed the parameters for 500 random projections of 
each galaxy and show the area enclosing the 90% most 
likely observed data points. For comparison we show the 
observations for disky (diamonds) and boxy (squares) el¬ 
lipticals galaxies. All galaxies are similar to isotropic 
rotators (defined here as log(v* ia j, ) > 0.7 ). Galaxy 
A is only marginally disky, Galaxy C can be very disky 
but also has some boxy projections, both are in good 
agreement with observations of disky and boxy isotropic 
ellipticals. Galaxy E has the largest values for 
and is extremely disky. It also show the largest values 
for rotational support v ma] / ao > 1 (Fig. [HI right panel) 
whereas galaxy C shows smaller values u ma j/<r 0 « 0.7 and 
ellipticites as large as e = 0.5. Galaxy a has the lowest 
values of v ma .i/a 0 ~ 0.2 and is nearly round ate = 0.2. 

In summary the high star-formation rate, the ex¬ 
ponential surface brightness distribution, the kine¬ 
matics and the isophotal shape of galaxy E makes 
it more similar to an early-type spiral galax¬ 
ies. Galaxies A and C do resemble IW* ellipti¬ 
cal galaxies. They follow the Fundamental Plane 
([ao,SB e ,r e ] = [153 krn/s, 20.5magarcsec -2 , 2.7kpc]; 

[op , SB P ., ry] = [145 km/s, 20.3magarcsec -2 , 2.8 kpc]) 

(| Bender et a.l JIT992I) and are consist ent with the observe d 
n — Mb and the n — a correlation dTruiillo et al .1 [20041 ). 

6 . DISCUSSION 



Radius [r/r e ] 

Fig. 13.— Sam as Fig. 1121 for galaxy C. Note that this galaxy 
has a counter-rotating core with a radius of s» 0.5r e . 



- 0.2 


- 2-10 1 2 
local v/ct 

Fig. 14.— Local correlation between /13 and vex/uo for galaxies 
A,C, and E for 500 radom projections. 90% of the data points 
cover the regions indicated by the shaded areas. The observed 
values for disky (open diamonds) an d boxy ellipticals (filled boxes) 
are overplotted IBender et al.l[T994h . 

We have followed the formation and evolution of three 
~ M* field galaxies using numerical simulations ’ab ini¬ 
tio’ from cosmological initial conditions at high redshift. 
All galaxies start to form when cold gas that has col¬ 
lapsed in subunits is effectively consumed by star for¬ 
mation. At early times the remaining diffuse gas is pre¬ 
dominantly heated by shocks. In general the stars galax¬ 
ies assemble both by in situ star formation as well as 
major/minor mergers and accretion. During the early 
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Fig. 15.— Left panel. Isophotal shape a4 e g versus anisotropy parameter (n ma j/o'o)*• The contours indicate the 90% probability to find 
a simulated galaxy (seen from 500 random viewing angles) in the enclosed area. The horizontal line at (u ma j/cro)* = 0.5 separates the 
anisotropic boxy from the more rotationally flattened disky and boxy ellipticals. Right panel: Ellipticity e e fj versus the ratio of major axis 
rotation velocity and central velocity dispersion {y ma j / ag). The data for observed disky (diamonds), very anisotropic boxy (filled boxes) 
and less anisotropic boxy (open boxes) ellipticals have been kindly provided by Ralf Bender. The properties of the galaxies are largely 
consistent with rotationally flattened disky and boxy ellipticals. 


formation phase at 2 < z < 8 the assembly of the galax¬ 
ies is dominated by mergers of gas rich subcomponents 
and in situ star formation. Although some fraction of 
the stars is accreted this phase has the characteristics 
of a dissipative collapse. Thereafter stellar accretion or 
minor/major mergers become more important and , for 
early-type galaxies, tends to dominate for z < 1. 

After z=l the further assembly of the galaxies varies 
from halo to halo. The evolution of one galaxy is domi¬ 
nated by gas infall and in situ star formation at rates of 
5 — 8 M©/yr. It develops a massive, however too compact, 
stellar disk and has properties more similar to early-type 
disk galaxies. 

Two other galaxies have present day properties more 
similar to elliptical galaxies and their assembly is dom¬ 
inated by minor and major predominantly stellar merg¬ 
ers. One of the galaxies grows by minor mergers and 
does not experience a major merger after its formation 
phase. At redshi ft z k, 1 it has formed a massive stellar 
spheroid (see e.g. iNagamine et all2005l ) which resembles 
observed massive evo l ved gala xies a nd their inferred for¬ 
mation histories dDaddi et al.ll2000f ). The system evolves 
into an old stellar system surrounded by a halo of hot 
gas at the present day. The second early-type galaxy 
had a major merger at z= 0.6 that can be classified as a 
gas poor or ’dry 1 merger as it was not accompanied by a 
strong burst of star-formation. From those two examples 
we find that ss 15% - 35% of the final stellar mass was 
assembled by accreti on of stars that have f ormed out¬ 
side the galaxy fsee iKhochfar fe Silkl l2006bl) . This in¬ 
crease will contribute to the evolution of the luminosity 
function of elliptical galaxies from z = 1 — > 0 in addi¬ 
tion to a fading population of blue galaxies and merg¬ 
ers (Bell et alJl2004 iDrorv et al.ll2004 IConselice et al.1 
12005 : iFaber et all 120051) . The size evolution of the sim¬ 
ulated early-type galaxies is determined by the accreted 
stars. While the half-mass radius of the stars formed in- 
situ remains almost constant independent of redshift at 
r i /2 = 1 — 2 kpc the accreted stars create an envelope 
whose hal f-mass radius is increasin g with decreasing red¬ 
shift (see IKhochfar fe Silk! l2006al for a semi-analytical 


approach). Interestingly, a significant fraction of the 
stars at the final half-mass radius of the galaxies were 
not born there but accreted. This process might pro¬ 
vide a good explan ation for a possible size evolution of 
early -type galaxies (iTruiillo et al.ll20(56t iLonghetti et al.l 
2006) without changing the fundamental parameter re¬ 
lations which are determined by the central properties of 
the galaxies and should have been established during the 
early formation phase. 

At the present day most investi gate d k inematical an d 
photometric properties (see also iSommer-Larsen et al.1 
120031) are consistent with observed rotating disky and 
boxy early-type galaxies at similar luminosities. In par¬ 
ticular, the concentration is in good agreement with 
observation s. W e believe that the discrepancy with 
iMeza et ahl (l2003f) is due to the different star forma¬ 
tion algorithm which in their case allows predominantly 
gaseous subunits to merge and develop a high concen¬ 
tration before forming stars. In our simulations gas is 
converted into stars efficiently in subunits which merge 
at early times, prior to z= 2 , and so the morphology of the 
systems is elliptical if its late accretion is dominated by 
stars. Interestingly, the early-type galaxy with the late 
major merger develops a counter rotating core which pre¬ 
viously only have been successfull y simulated in isolate d 
equal-mass disk-di sk mergers dHernauist fe Barnesll 19911 
iJesseit et al.l 120061) and are supposed to be signposts of 
merger events. 

Feedback from supernovae and/or AGN, which was 
not included in this simulation, would provide additional 
heat sources but appear neither vital to stop the star 
formation and produce a massive red galaxy at higher 
redshifts nor to produce a present day stellar system 
with consistent kinematical properties. Theoretical argu- 
ments (see e.g. iDekel fe Silklll986t lEfstathioull2000D sug¬ 
gest that stellar feedback rapidly becomes ineffective in 
systems with velocity dispersions greater than 100 km/s. 
However, AGN feedback may be important in galaxies 
with high velocity dispersion s (|Silk fe Reeslll998t iFabianl 
ll999tlCiotti fe Qstrikeill200ll) . In our case AGN feedback 
might have been able to prevent further star formation 
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at low redshift by heating the cold gas at the center of 
the g alaxy (see e.g. ISoringel et all[200 5 al ISazonov et al.l 
2005). On the other hand many field ellipticals show the 
related E+A phenomenon which seems to be suppressed 
in the centers of clusters as any re sidual cold gas will 
be blown away ( Caldwell et al.lll993h . For more massive 
systems than the ones studied here it has been argued 
that feedback from AGN might have a bigger impact and 
be responsible for the sharp cutoff at the bright end of 
the luminosity function of ellipticals (|Croton et al .Il2005h . 
However, only full cosmological simulations at high res¬ 
olution of the formation of more massive systems will 
help to fully understand the importance of stellar and 
AGN feedback on the formation and evolution of ellipti¬ 
cal galaxies. 

In a hierarchical universe most galaxies are expected 
to experience mergers but it has been argued for a 
long time that dissipation is essential to explain prop¬ 
erties of elliptical galaxi es an d _in partic ular inter¬ 
mediate mass elli pticals (Kormend y 19891 and refer¬ 
ences therein; iKormendv fe Bender! 19961 : iNaab et alJ 
l2006al) . The cosmological formation processes pre¬ 


sented here combines characteristics of both stellar merg¬ 
ers and dissipational collapse/mergers whereas the lat¬ 
ter dominates the early formation phase. It natu¬ 
rally meets many observational constraints on the evo¬ 
lutionof massive spheroidal galaxies (e.g. iTreu et al.l 
120051 : IThomas et aLll2005at Ivan der Wei et, al.ll2005l l. It 
might be equally or even mor e important than ma¬ 
jor mergers of dis k galaxies (R obertson et al.l 120061 : 
iHookins et al.ll2006l ) which can n ot reproduce the prop¬ 
erties the most m assive ellipticals (INaab fe Burkertll2003t 
I Cox et all I2006H whose evol ution could be driven by 
mergers of evolved ellipt icals (Ostrikcr & HausmaiJlHIi 
Khochfar fc BurkertllAiO.'ll iBell et al. 2006a; Naab et all 
2006briBovlan-Kolchin et al.ll2006lf 
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able prior to publication and Ralf Bender and Andy 
Fabian for discussions about the manuscript.We also 
thank the anonymous referee for valuable comments on 
the manuscript. 
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